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Abstract 

We update the branching ratios for the inclusive decays B —> and the 

exclusive decays B —> {K, , with £ = e, /r, in the standard model by including 

the explicit 0{as) and Agcn/m-fc corrections. This framework is used in conjunction 
with the current measurements of the branching ratios for B —> Xgj and B KBi~ 
decays and upper limits on the branching ratios for the decays B —> (K*, Xs)B£~ to 
work out bounds on the Wilson coefficients Cj, Cg, Cg and Gio appearing in the effective 
Hamiltonian formalism. The resulting bounds are found to be consistent with the 
predictions of the standard model and some variants of supersymmetric theories. We 
illustrate the constraints on supersymmetric parameters that the current data on rare B 
decays implies in the context of minimal flavour violating model and in more general 
scenarios admitting additional flavour changing mechanisms. Precise measurements 
of the dilepton invariant mass distributions in the decays B —> {Xg, K*, K)B£~, in 
particular in the lower dilepton mass region, and the forward-backward asymmetry in 
the decays B —> {Xg, K*)B£~, will greatly help in discriminating among the SM and 
various supersymmetric theories. 
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1 Introduction 


The measurement of the inclusive decay B —>■ Xg'j, hrst reported by the CLEO collaboration 
in 1995 [Q], has received a resounding reception in the interested theoretical physics commu¬ 
nity, both as a precision test of the standard model in the flavour sector and as a harbinger 
of new physics, in particular supersymmetry |^. In the meanwhile, the branching ratio for 
this decay has become quite precise through the subsequent measurements by the CLEO , 
ALEPH m and BELLE collaborations, with the BABAR measurements keenly awaited. 
The present experimental average of the branching ratio B{B —> = (3.22 ±0.40) x 10“^ 

is in good agreement with the next-to-leading order prediction of the same in the standard 
model (SM), estimated as B{B —> Xs7)sm = (3.35 ± 0.30) x 10““^ |^,|^ for the pole quark 
mass ratio rric/mi, = 0.29 ± 0.02, rising to B{B —>• Xs7)sm = (3.73 ± 0.30) x 10““^ |@], if 
one uses the input value rri^^{fi)/rrib^poie = 0.22 ± 0.04, where m^®(/i) is the charm quark 
mass in the MS-scheme, evaluated at a scale /i in the range rric < fi < mb. The inherent 
uncertainty reflects in part the present accuracy of the theoretical branching ratio, which is 
limited to 0{as), and in part the imprecise measurements of the photon energy spectrum 
in R —>■ Xg'y decays. Despite the current theoretical dispersion on the branching ratio, the 
agreement between experiment and the SM is quite impressive and this has been used to put 
non-trivial constraints on the parameters of models incorporating beyond-the-SM physics, 
in particular supersymmetry (see, for example. Ref. |P] for a recent analysis in a supersym¬ 
metric scenario). While the measurement of the decay B X^y is being consolidated, 
several other radiative and semileptonic rare R-decays are being searched for. In partic¬ 
ular, first measurements of semileptonic rare R-decays have been recently reported in the 
R —> and R — Ke^e~ modes by the BELLE collaboration and upper limits 

have been put in a number of other related decay modes |l^-[l^. The current measurements 
of the exclusive modes are in agreement with the expectations in the SM |T^,0, calculated 


in next-to-leading logarithmic (NLO) approximation, taking into account the experimental 
and theoretical errors. This, for example, can be judged from the comparison of the com¬ 
bined branching ratio for the decay modes R ^ K£^£~, £ = e,fi, reported by the BELLE 
collaboration B{B —>■ K£^£~) = 0.75lo;2i ± 0.09 x 10“® with the Light-cone QCD sum rule 
based estimates of the same, B{B —^ K£^£~) = 0.57to;}o x 10“® 
inclusive decays R ^ Xs£'^£~ and the exclusive decays R —> K*£~^£ 


The upper limits on the 
are now approaching 


their respective SM-based estimates, as we also show quantitatively in this paper. 

With increased statistical power of experiments at the R factories in the next several 
years, the decays discussed above and related rare R decays will be measured very precisely. 
On the theoretical side, partial results in next-to-next-to-leading logarithmic (NNLO) accu¬ 
racy are now available in the inclusive decays R —> Xs£~^£~ - BOI- Recalling that the lowest 
order contribution for these decays starts at 0{l/as), as opposed to the decay R —>• X^y, 
which starts at 0{a^), the NNLO accuracy in R — Xs£~^£~ amounts to calculating explicit 
0{as) improvements. The same accuracy in as amounts to calculating the decay R —X^y 
in NLO. We also recall that power corrections in Aqqd/ mb []^ and in Aqcd/'^c |0 are 
also known. What concerns the exclusive decays, some theoretical progress in calculating 
their decay rates to NLO accuracy in the R —> (X*,p)y , and to NNLO accuracy in 
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B —> [^] decays, including the leading Aqcb/Mb, has been reported. Comparisons 

of these theoretical estimates with data on B ^ K*'y decays have led to important 


inferences on the magnetic moment form factor. Our purpose in this paper is to incorporate 
these theoretical improvements, carried out in the context of the SM, and phenomenologi¬ 
cal implications from the observed radiative decays, and examine the quantitative rapport 
between the SM and current measurements of the semileptonic rare 5-decays. An equally 
important undertaking of our analysis is to investigate the impact of the current experimen¬ 
tal measurements on the parameters of the possible supersymmetric extensions of the SM. 
The question which we address in this context is the following: do the current measurements 
in semileptonic rare 5-decays already provide more restrictive constraints on the parameters 
of the supersymmetric models than are provided by the 5 —> measurements? We hnd 

that the decays 5 —>■ (X^, X*, do provide additional constraints in some parts of the 

supersymmetric space, though with the current experimental knowledge the decay 5 —»■ X^q 
remains more restrictive over most of the supersymmetric space. This is expected to change 
with improved precision on the semileptonic rare 5-decays, which we illustrate in a number 
of supersymmetric scenarios. 


Our analysis is carried out in the effective Hamiltonian approach, obtained by integrat¬ 
ing out the heavy degrees of freedom, dehned below (see Eq. (|^). However, we make the 
tacit assumption that the dominant effects of an underlying supersymmetric theory can be 
implemented by using the SM operator basis for the effective Hamiltonian. Thus, super- 
symmetric effects enter in our analysis through the modihcations of the Wilson coefficients 
which in the SM are calculated at some high scale, denoted generically by piy, with the 
SM anomalous dimension matrix controlling the renormalization of these coefficients to a 
lower scale, typically Hb = 0{mb). Restricting the operator basis to the one in the SM 
obviously does not cover the most general supersymmetric case, but we think that it covers 
an important part of the underlying parameter space, and hence can be employed to un¬ 
dertake searches for supersymmetric effects in rare 5-decays. Within this operator basis, 
we have split our analysis in two parts. In the hrst part we update the branching ratios 
for the decays 5 — {Xg, K, , i = e, /i, in the standard model. In doing this, 

we work out the parametric uncertainties due to the scale-dependence, top quark mass, 
rrit, and the ratio of the quark masses rric/rab. Combining the individual errors 55(/i), 
5B{mt), and 6B{mc/mb) in quadrature, we hnd that the resulting theoretical uncertainties 
are 6B{B —>• XsC+e”) ~ ±15% and 6B{B —>• Xs/i+/r“) ~ ±17%. The corresponding theoret¬ 
ical uncertainties on the exclusive decay branching ratios are larger, due to the form factors, 
and estimated at typically 0(±35%). Using this updated theoretical framework, we extract 
mo del-independent constraints that current data (summarized below in Eqs. (|^)-(P%|)) pro¬ 
vides on the Wilson coefficients Cj - Cio, which appear in the effective Hamiltonian. We 
hrst work out the constraints on Cj and Cg implied by the 5 —X^q measurement. To that 
end, we dehne the quantities RY^siAw) = /Ci, and work out bounds on them. 

Data on 5 —>• X^q allows both Rj^siAw) > 0 and < 0 solutions, which we show in 

terms of the allowed regions in the {Rj{piw),R?,{.Aw)) and {Rj{pib),R&{Ah)) planes. We then 
transcribe the impact of the 5 —> (X^, K, K*)£^£~ experimental data on the allowed regions 
in the [Cg, Cio] plane. Depending on the two branches for the quantities Rj^si^w), we display 
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the constraints in terms of Cg^{fiw) and We show that the SM solntion (corresponding 
to the point (0, 0) in this plane for the case = 1) is allowed by present data. More 

importantly, from the point of view of snpersymmetry, onr analysis shows the allowed region 
in the [Cg, Cio] plane, which leaves considerable room for beyond-the-SM contribntions to 
these qnantities. In fact, in some allowed regions, phenomenological prohles of semileptonic 
rare i?-decays can measnrably differ from the corresponding ones in the SM. 


The second part of onr snpersymmetric analysis deals with specihc SUSY models and 
we qnantify the additional constraints that the generic b —>■ si^i~ data implies for the 
parameters of these models. We show that no nsefnl bonnds beyond what are already known 
from the B —> Xgj analysis are at present obtained in the so-called minimal flavonr violating 
models (inclnding the constrained minimal snpersymmetric standard model MSSM) . 
This reflects the generically small deviations to the SM rates and distribntions anticipated 
in these models, as the allowed snpersymmetric parameters are already highly constrained. 
Working in the mass insertion approximation |^, we show that insertions in the down- 


sqnark sector (that enter principally throngh the glnino pengnin and box diagrams) are not 
constrained either from present data. On the other hand, insertions in the np-sqnark sector 
get in some parts of the SUSY parameter space genninely new constraints. To show possible 
snpersymmetric effects that precise measnrements in semileptonic rare 5-decays may reveal, 
we work ont the forward-backward asymmetry in 5 —>■ for fonr illnstrative points in 

the ( 0 |^^(pvv )5 plane, representing solntions in the fonr allowed qnadrants in this space. 
However, a high density scan over all the parameter space shows that the allowed solntions in 
the models considered by ns are scattered mostly aronnd the (0|^^(/ivv), = (0, 0) region, 

for the two branches for the qnantities i? 7 , 8 (hw), be. for > 0 and RrfiifJ'w) < 0. 

We present the resnlting constraints on the snpersymmetric masses (mass of the lighter 
of the two stop mass eigenstate), Mh± (the charged Higgs boson masses) and tan/? (ratio 
of the two Higgs vacnnm expectation valnes), in the context of the MFV-MSSM framework, 
and on the mass insertion parameter (^ 23 ) in the MIA framework. This npdates similar 
resnlts worked ont along these lines in Ref. [^ . 

This paper is organized as follows: In Section we enlist the cnrrent measnrements of 
the rare 5-decays which we have analyzed. The effective Hamiltonian for the SM and the 
snpersymmetric models stndied by ns is also given here. In Section |^, we present the NNLO 
implementation of the inclnsive and exclnsive b —»• transitions that we consider. In 

Section ^ we discnss the branching ratios for the exclusive decays 5 —in the SM. 
In Section we study the constraints on the supersymmetric parameters resulting from the 
5 —> Ys 7 decay in the NLO approximation. In Section we present the results of the 
model-independent analysis of the b —>• st^l~ modes based on current data. In Section |^, 
we describe the specihc SUSY model that we study and present the bounds on the relevant 
mass insertions. In Section we summarize our results. Some loop-functions encountered 
in the calculation and the stop and chargino mass matrices are given in the appendices. 
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2 Effective Hamiltonian 


The effective Hamiltonian in the SM inducing the b —»• si~^i and b ^ sj transitions can be 
expressed as follows: 


^eff 


4Gf 

■ 71 “ 


10 


K:VaY.C,{p‘)0.{p) 


( 1 ) 


2=1 


where Oi{fi) are dimension-six operators at the scale /r, Cj(/i) are the corresponding Wilson 
coefficients, Gf is the Fermi coupling constant, and the CKM dependence has been made 
explicit. The operators can be chosen as Ref. 


16 


Oi 

O 3 
O, = 


= {sL'y^T^CL){cLYT'^h), 

= {sL7tibL)Eq{qYq), 


O 2 

O 4 


{sL7^,CL){cLYh) , 

(sL^t^T^bL) EgiqYT'^q) , 

(sL7Mi7M2 7/.3^“^i) Eq{qY"Y^Y^T<^q) , 


O7 = ■^mb{sL(J^^’'bR)F^ 


flU ? 


Os = ^mb{sL(T^^^T'^bji)Gl 


O 9 = 


Oio — 


^isL7^^bL)Eei^7^75i) , 


( 2 ) 

where the subscripts L and R refer to left- and right- handed components of the fermion helds. 
We work in the approximation where the combination (IG Wb) of the Cabibbo-Kobayashi- 


Maskawa (CKM) matrix elements [^] is neglected; in this case the CKM structure factorizes, 
as indicated in Eq. (|^. Of course, for the sake of book keeping, one can keep the individual 
top-quark and charm-quark contributions in the loop separately, but as there is no way to 
distinguish these individual contributions we will give the results in the summed form. 

Note the inverse powers of gs in the definition of the operators 07 ,..., 0 io in Eq. (2). These 
factors have been introduced by Misiak in Ref. |^] in order to simplify the organization of 
the calculations. In this framework, the LO result for the b —> s£~^i~ decay amplitude is 
obtained in the following three steps: the matching conditions Gi{fiw) have to be worked 
out at 0(q:°), the renormalization group evolution has to be performed using the 0{al) 
anomalous dimension matrix and the matrix elements of the operators Oj have to be worked 
out at order l/a<j. In going to the NLO precision all the three steps have to be improved by 
one order in ag- 

At an arbitrary scale /i the Wilson coefficients can be decomposed as 


c,{p) = cf (ri + 


asiqi^) 

Att 


G, 




ifi) + 


(3) 


We note that in our basis only G 2 is different from zero at the matching scale pw at leading 
order, viz. G^^\iJ,w) = ^i 2 - At the low scale pb (of order mb), the coefficients G^{iJ,b) are 
non-zero for i = 1,..., 6; 9 whereas they vanish for i = 7,8, 10. 

We shall use this effective Hamiltonian and calculate the matrix elements for the decays 
of interest, specifying the degree of theoretical accuracy. 
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The experimental input that we use in our analysis is given below. Except for the inclusive 
branching ratio for B —*■ which is the average of the results from CLEO, ALEPH and 

BELLE measurements p-|^, all other entries are taken from the two BELLE papers listed 
in Ref. [|T^ : 


= (3.22 ± 0.40) X 10"^ , (4) 

= x 10“^ , (5) 

= X 10 -^ , ( 6 ) 

= ± 0.09) x 10"® , (7) 

< 3.0 X 10"® at 90% C.L. , ( 8 ) 

< 5.1 X 10"® at 90% C.L. , (9) 

< 19.1 X 10“® at 90% C.L. , (10) 

< 10.1 X 10“® at 90% C.L. . (11) 

IS. ® - 0 ) refer to the so-called non-resonant 
branching ratios integrated over the entire dilepton invariant mass spectrum. In the ex¬ 
perimental analyses, judicious cuts are used to remove the dominant resonant contributions 
arising from the decays B —> (X^, K, ,...) —>• (X^, K, A direct compar¬ 

ison of experiment and theory is, of course, very desirable, but we do not have access to this 
restricted experimental information. Instead, we compare the theoretical predictions with 
data which has been corrected for the experimental acceptance using SM-based theoretical 
distributions from Ref. |T^,|^ . In the present analysis, we are assuming that the acceptance 
corrections have been adequately incorporated in the experimental analysis in providing the 
branching ratios and upper limits listed above. We will give the theoretical branching ratios 
integrated over all dilepton invariant masses to compare with these numbers. However, for 
future analyses, we emphasize the dilepton invariant mass distribution in the low-s region, 
s = 771^+^-/ml < 0.25, where the NNLO calculations for the inclusive decays are known, 
and resonant effects due to J/i/, 'i// etc. are expected to be small. 


B{B ^ X,7) 
B{B Xp+/i-) 
B{B Xe+e") 
B{B Ktr) 
B{B XVh“) 
B{B K*e+e-) 
B{B ^ X,/r+/i-) 
B{B X,e+e-) 

The experimental numbers given in 1 


3 Inclusive b si^i decays at NNLO 


We start by discussing the NNLO analysis of the B —> Xgi^i decays presented in Refs. 


T7| , recalling that the 0(a<j) corrections to the matrix elements computed in Ref. [|^ have 
been calculated only below the cc resonances. 

In the NNLO approximation, the invariant dilepton mass distribution for the inclusive 
decay B —> can be written as 


dTjb ^ x/+r) _ 

ds \ 47r 

(1 -|- 2s) ( Ct 


2 n2 5 


eff 


+ 


C- 


eff 

10 


4871 ^ 

+ 4(1 + 2/s) c: 


—(1 - s) 


X 


eff 


+ 12Re ( 


( 12 ) 
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In the SM the effective Wilson coefficients C'®®, and are given by 


^eff 


= 


= 


/^eff _ 
*-^10 ~ 


, «.(h) /.A . 

1 H-INTI'S) I Ay 


TT 


^ (C^Ff)(S) + Ff (s) + 4°^Fi")(s)) , 

1 H- ^^UJg{s)^ (^Ag + Tg h^TTl^, s) + Ug h(l, s) + Wg h(0, s)^ 

^ (F^F® (s) + Ff Ff (s) + A^i'^F^^\s)) , 


4:71 

'l + ^^s(s)] ^10 , 


(13) 

(14) 

(15) 


where the functions h{ml,s) and a; 9 (s) are given in Ref. |^, while ci; 7 (s) and F^^ 2 , 8 H'S) can 
be seen in Ref. [^. The auxiliary quantities R 7 , A^, Ag, Aig, Tg, Ug, ITg are the following 
linear combinations of the Wilson coefficients Ci{fi) (see Eq.(|l|)): 


-47 = ^C 7 (,i)-lc 3 (i.)-ic 4 (i.)-|!c 5 (,.)-|!c«M, 

.^8 — -7—7 Csip) + C^{p.') — — Ci{p.') + 20 C,{^) —— Cg{^) , 

as[^i) 0 6 

^9 = Cg{fi) + ^ Cii^i) 7 ^ 9 ^ In — 

Ots\IU jJ, 


(16) 

(17) 



4 64 64 

+ 2 U3(/i) + — F5(/i) -h — Ceifi), 

(18) 

^10 

477 

(M) 

(19) 

Tg 

= -h--h 02 (^ 7 ) -h 6 -h 60 Cslfl) , 

( 20 ) 

Ug 

7 2 32 

= --F3(/x)--F4(/i)-38F5(/i)-yF6(/i), 

( 21 ) 

Wg 

12 32 

= --F3(/i)--F4(/i)-8C'5(/i)-yC'6(/i). 

( 22 ) 


The elements j-g^ can be seen in Eq. (26) of ref. in Eqs. (p!3[) and (|T^ denotes the 

lowest order piece of Ag: 


4“’ = d‘’(fi) + cf (a) - icfF) + 20 Gf (fi) 

6 

(0) 


10 ^(0) 


cr(M) • 


(23) 


The numerical values for the coefficients 247 ,^ 3 ,^ 49 , Aig, Tg, Ug, Wg, Ci and Cg are obtained 
after solving the renormalization group equations for the Wilson coefficients Ci{fi), using the 


matching conditions from ref. and the anomalous dimension matrices from Refs. 

As mentioned earlier, we do not separate charm- and top- quark contributions and perform 
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the matching (for both) at the scale /ivy = 'n^w- The resulting values are summarized in 
Table 01. Note that when calculating the decay width (IT^), we retain only terms linear in 
a, (and thus in oig and uj^) in and In the interference term Re 

too, we keep only terms linear in By construction, one has to make the replacements 
Ug —>• cuyg and ujy —>• a ;79 in this term where the function cj 7 g(s) can be found in Ref. . 

We now turn to the modifications of the effective Wilson coefficients , Cg and Cfg 
in the presence of new physics which enters through a modification of the Wilson coefficients 
Cy, Cg, Cg and Cio at the matching scale By doing so, we tacitly assume that the scale of 
new physics is close enough to the weak scale mw, justifying to integrate out simultaneously 
the heavy SM particles and the additional ones present in the new physics scenario. For 
simplicity we assume that only the lowest non-trivial order of these Wilson coefficients get 
modified by new physics, which in our setup (see Eqs. ®,i) ,(11)) means that Cj^^jUw), 
Cg^\jUiv), Cg^^jUiv), CiQ^jUiy) get modified. The shifts of the Wilson coefficients at jUiy can 
be written as: 

Ciifiw) —^ ■ (24) 

47r 

These shift at the matching scale are translated through the RGE step into modifications of 
the coefficients Ci{fih) at the low scale leading in turn to modifications of the effective 
Wilson coefficients defined in Eqs. (pT^-jlSD. They now read 


cf = 


F'cff _ 

Ug — 


r^eff 


as{ii) 


1 H- - — oj7(s)j (Ay + Ayy Cy (i^w) + ^78 Cg (l^w)) 

(jC) 


(cfFP{g) + cf + 4">Ff>(S) + 4? FP{g)) ,(25) 

(yl, + T, j) + u, ;,(1, j) + Wg h(0, s) + C^'"(pw)) 

^ (Ci”'q’>(s) + cf Ff >(S) + Af q’>(S) + 4”' Crimr) d”(s)) .(26) 

( 27 ) 



1+^^9(i)) (21,„ +csn. 


The numerical values for the parameters Ryy, Ayg, Rgg^ which incorporate the effects from 
the running, are listed in Table |^. 


3.1 Power corrections in inclusive B 


X,t^t 


decays 


Before presenting a theoretical analysis of the available data on rare R-decays, we would like 
to discuss power corrections in the inclusive B —> XsBl~ decays. In the NNLO approxima¬ 
tion and including leading order power corrections in l/rub [Q and l/rUc 
dilepton mass distribution for the inclusive decay B XsBl~ can be written as 


191, the invariant 


dT{b si+t 
ds 


OLp 


Itt 


2 G\ 


m 


b,pole 




tb I 


4877^ 


■(1-5)^ 


(1 -f 2s) 


G 


eff 


+ 


C 


eff 

10 


Gi(s) 
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fi = 2.5 GeV 

/i = 5 GeV 

fi = 10 GeV 

(y^ 

0.267 

0.215 

0.180 

(cr.cp) 

(-0.697,0.241) 

(-0.487,0.207) 

(-0.326,0.184) 

X 3 

o' 

(1.046,-0.028) 

(1.024, -0.017) 

( 1 . 011 ,- 0 . 010 ) 

(7ir,<') 

(-0.353,0.023) 

(-0.312,0.008) 

(-0.278,-0.002) 

Mxwr 

(0.577, -0.0524) 

(0.672,-0.0391) 

(0.760, -0.0277) 

(aW 4 + ^ 
12^78 , /I 78 J 

(0.109,-0.00520) 

(0.0914, -0.00193) 

(0.0707, -0.000263) 


-0.164 

-0.148 

-0.134 

2^88 

0.618 

0.706 

0.786 

(+’.+’) 

(4.287,-0.218) 

(4.174, -0.035) 

(4.177,0.107) 

■p+T+j— 

(0.114,0.280) 

(0.374,0.252) 

(0.575,0.231) 

(+’.+’) 

(0.045,0.023) 

(0.033,0.015) 

( 0 . 022 , 0 . 010 ) 

IWXfS 

(0.044,0.016) 

(0.032,0.012) 

(0.022,0.008) 

(aW p + n 
1^10 ) 2^10 j 

(-4.592,0.379) 

(-4.592,0.379) 

(-4.592,0.379) 


Table 1 : Coefficients appearing in Eqs. (|T3|-|T^) and Eqs. for three different scales 

fi = 2.5 GeV, fi = 5 GeV and /r = 10 GeV. For as(p) (in the MS scheme) we used the 
two-loop expression with 5 flavours and as{mz) = 0.119. The entries correspond to the MS 
top quark mass renormalized at the scale mw, mt{mw) = 175.9 GeV. The superscript (0) 
refers to the lowest order quantities while the superscript ( 1 ) denotes the correction terms 
of order t.e. X = with X = C, A, T, U, W. 


where 


4 (1 + 2/s) 


a 


eff 


G2(3) + 12Re ( Cfcf *) ^ 3 ( 5 ) + G,(s) 


— 1 _L 1 — 15s^ + lOs^ A 2 

Gi(s) — 1 + 7;—7 + 3 


2ml (1 — s) 2 (l + 2 s) 2 m^ 
3 


('■\ 1 I 6 + 3s — 5s A 2 

G 2 s = 1 + 7 ^-3 


G3(s) = 1 + 


2ml (1 — s)2(2 + s) 2ml 
Ai 5 + 6s — 7s^ A 2 


2ml 


(1 


2m 


2 • 


(28) 

(29) 

(30) 

(31) 


The values of the heavy quark matrix elements Ai and A 2 that we use in our analysis are 
given in Table 2. The term denoted by Gc takes l/rric corrections into account. It is written 
as 






(32) 


Since our basis in Eq. (^) is different from the one often used in the literature i.e. Oi = 
('SL 7 /.t&L)(cL 7 ^Ci) and Oi = Oi/2 —O 2 / 6 , Eq. (|3^) differs superhcially from the one reported 



































in The function F{r), where r = s/(4m^), is given below ||19|| : 

1 


F(r) = — < 
^ ^ 2 r ' 




arctan. 


1 — r 


- 1 


1 - Jl - 1 /r 

In- ^, + ivr I — 1 


, 2^r{r - 1) V 1 + - 1/ 


0 < r < 1 , 


r > 1 . 


(33) 


The impact of power corrections in inclusive decays B —> Xsi^i~ at NLO has been 
studied in [^] in the SM. In the low dilepton mass region and for not too close to 
the photon pole where O 7 dominates, the 1/mb effects enhance the rate by ~ 1%. In the 
high-s region they become negative and decrease the rate by few percent. Their magnitude 
rises more and more towards the boundary ~ mf, where the expansion in l/rub breaks 
down 1]^. The 1/mc expansion Eq. 


is valid everywhere except near threshold s = 4m^, 
and it also fails at the charmonium resonances J/T and higher ones like 4/'. The 1/mc 
corrections decrease the rate below the charm threshold and enhance it above by few percent. 

This is illustrated in Figure ^ where the relative size R{s) of the combined l/rrib and 
1 /mc corrections dehned as 


m = 


^ (with power corrections) — power corrections) 


dr{b- 


■'Xsi+e-) 

ds 


(with power corrections) 


(34) 


is shown for the SM, and for comparison also for C^ = —Cj^. Both 1/m correction thus 
partially cancel in the SM. The situation with new physics can be different. In a generic 
scenario with C 7 = —Cj^ the power corrections can be more pronounced, in particular for 
low dilepton mass where both 1/m corrections are negative. Together they lower the rate 
by few percent. Note that in our estimates of the B —>■ branching ratio, we include 

the power corrections in the semileptonic branching ratios [0 as well. 


3.2 Branching ratios for B Xs£~^£ in the SM 


In order to eliminate the large uncertainty due to the factor appearing in the decay 

width for B —>■ XsBi~, it has become customary to consider instead the following branching 
ratio 




dT{B XJ+j-) 
r(S —>■ XcCp) ds 


(35) 


in which the factor drops out. The explicit expression for the semi-leptonic decay 

width r(i? — XcCz/g) can be found e.g. in Ref. |^. Note that as we are ignoring the 
annihilation contributions, which lead to isospin violations in the decay widths, and we are 
using the averaged semileptonic branching ratio to normalize, all our inclusive branching 
ratios are to be understood as averaged over the B^ and B^{B^) decays. 


The dilepton invariant mass distribution for the process B —> calculated in 

NNLO is shown in Fig. ^ for the three choices of the scale /r = 2.5 GeV, /i = 5 GeV and 
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fx = 10 GeV (solid curves). In this figure, the left-hand plot shows the distribution in 

the very low invariant mass region (s G [0,0.05], with 0 to be understood as the kinematic 
threshold s = drUg ~ 10“® GeV^, yielding s = 4.5 x 10“®), and the right-hand plot shows 
the dilepton spectrum in the region beyond s > 0.05, and hence this also holds for the 
decay B —> Xsp+/i“. We should stress at this point that a genuine NNLO calculation only 
exists for values of s below 0.25, which is indicated in the right-hand plot by the vertical 
dotted line. For higher values of s, an estimate of the NNLO result is obtained by an 
extrapolation procedure discussed in more detail at the end of this paragraph. The so-called 
partial NNLO dilepton spectrum, obtained by switching off the quantities ^ 1^23 Eqs. (^) 
and (|T^), is also shown in each of these cases for the same three choices of the scale /i (dashed 
curves). Note that in the left-hand plot, the lowest lying curves are for /i = 10 GeV and 
the uppermost ones are for fi = 2.5 GeV. In the right-hand plot, the scale-dependence is 
reversed, i.e., the highest lying curves are for p = 10 GeV and the lowest for /i = 2.5 GeV. 
The crossing (in the partial NNLO BR) happens near s = 0.04 and this feature leads to 
a certain cancellation of the p dependence in the decay rate for B —> Xse'^e~. We also 
note that the NNLO dilepton invariant mass spectrum in the right-hand plot (s > 0.05) 
lies below its partial NNLO counterpart, and hence the partial branching ratios for both 
the B —> Xse^e~ and B Xs/r+/i“ decays are reduced in the full NNLO accuracy. More 
importantly, from the point of view of our subsequent analysis. Fig. shows that the full 
NNLO invariant mass distribution is very well approximated by the partial NNLO for the 
choice of the scale p = 2.5 GeV, in the entire low-s range. This is yet another illustration 
of the situation often met in perturbation theory that a judicious choice of the scale reduces 
the higher order corrections. From this observation, it seems reasonable to use the partial 
NNLO curve corresponding to ph = 2.5 GeV as an estimate for the central value of the full 
NNLO for s > 0.25. We estimate the scale dependence in this region by assuming that it is 
given by the genuine NNLO calculation at s = 0.25. 

In order to complete our discussion of the computation of the inclusive branching ratios, 
it is necessary to discuss their dependence on the quark masses rrit and rric (in particular, 
the latter is marked for what concerns the SM prediction). We vary both masses within the 
errors that we quote in Table ^ and present the results for the branching ratios B —>■ Xse'^e~ 
and B W/x+/i“ in Table |] where we include also the power corrections discussed in 
Sec. |3.1| . In Table ^ we show the SM central values and the parametric uncertainties by 
means of independent error bars (to be interpreted as 68% G.L. uncertainties). In this table, 
the first error on the exclusive channel is due to the form factors, and is by far the dominant 
one. The other errors in both the exclusive and inclusive decays come from the scale (pfe), 
i^t,poie and mc/rrib respectively. Summing the errors in quadrature we get for the inclusive 
decays: 

B{B X,e+e“) = (6.89 ± 1.01) x 10“® {5Bxsee = ±15%) , (36) 

B{B ^ X,p+/i-) = (4.15 ± 0.70) X 10“® {SBx^pp = ±17%) . (37) 

Using the same input parameters, but restricting to the NLO precision, the inclusive branch¬ 
ing ratios have the central values B{B —> Xge+e") = 7.8 x 10“® and B{B —>• Xsfi~^fi~) = 
5.2 X 10“®. Thus, NNLO corrections reduce the branching ratios by typically 12% and 20%, 
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respectively. In Ref. 0, it was recently suggested in the context of the decay B Xg'j, 
where the charm quark mass enters the matrix elements at the two-loop level only, that it 
would be more appropriate to use the running charm mass evaluated at the pb ~ 0{mh) 
scale, leading to rric/mb — 0.22. Intuitively, this is a reasonable choice since the charm quark 
enters only as virtual particle running inside loops; formally, on the other hand, it is also 
clear that the difference between the results obtained by interpreting rric as the pole mass or 
the running mass is a higher order QCD effect. In what concerns B —>■ , the situation 

is somewhat different, as the charm quark mass enters in this case also in the one-loop matrix 
elements associated with Oi and 02- In these one-loop contributions, rric has the meaning 
of the pole mass when using the expressions derived in Ref. 
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Concerning the charm 

quark mass in the two-loop expressions, the dehnition rric is not hxed, like in R —>■ In 

our analysis, we prefer not to include this effect related to the dehnition of the charm quark 
mass in the hnal errors that we have listed. 


4 Exclusive B decays 


For what concerns the exclusive decays B —^ we implement the NNLO corrections 

calculated by Bobeth et al. in Ref. [|T^ and by Asatrian et ah in Ref. |]^ for the short-distance 
contribution. Then, we use the form factors calculated with the help of the QCD sum rules 
in Ref. |]^. Note that, in this case, we have dropped the contribution to the matrix elements 
given by the functions a;i(s) since this can be regarded as included in the full QCD form 
factors. In adopting this procedure, we are ignoring the so-called hard spectator corrections, 
calculated in the decays B —^ pl| in the large energy limit of QCD 1^, necessarily 


limiting the invariant mass to the small-s region. The hndings of Ref. [^| are that the 


dilepton invariant mass distribution in this region is rather stable against the explicit 0{as) 
corrections, and the theoretical uncertainties are dominated by the form factors and other 
non-perturbative parameters specihc to the large-energy factorization approach. This is 
so, even if one takes the point of view that the form factor ^±( 0 ), governing the transition 
B —> K* to the transversely polarized iF*-meson, can be assumed known from the analysis of 
the radiative transition B —>■ K*'^ in this approach and current data, as it is the contribution 
of the longitudinally polarized K* which dominates the decay rate in the small-s range, for 
which a knowledge of .^n is required. In principle, using estimates of SU(3)-breaking and 
HQET, the function ,^|| for the decays B can be obtained from the semileptonic 

decays B —>• piui. However, as present data on the Q^-dependence in the decay B —>• is 

rather sparse and a helicity-based analysis of the decays B —>■ ptv^ has yet to be undertaken, 
one will have to resort to form factor models for ,^||, which as opposed to the transverse form 
factor is essentially unbounded. In view of this, we ignore the hard spectator correction 
and discuss a plausible range of the form factors in the decays B {K, 


As already stated, some inference about the magnetic moment form factor Ti(0), involv¬ 
ing the matrix element of the operator O 7 in the decay B 77*7, has been derived by 
comparing the explicit 0{as) and Aqcd/Mb corrected branching ratios in the factorization 
approach with data pU|-[2^. One hnds that present data on R — 77*7 decay yields typically 
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a value in the range Ti(0) = 0.28 ± 0.04. This suggests that, including the explicit 0{as) 
corrections, data requires a value of this form factor which is smaller than its typical QCD 
sum rule estimate. To accommodate this, we use the minimum allowed form factors obtained 
in the light-cone QCD sum rule formalism, given in Table 5 of Ref. [|^, as our default set. 
This, for example, corresponds to setting Ti(0) = 0.33. In our numerical analysis, we add a 
flat ±15% error as residual uncertainty on the form factors. Thus, the input range for Ti(0) 
in our analysis T'i(O) = 0.33 ± 0.05 overlaps with the phenomenologically extracted value in 
the factorization approach given earlier. Again, following the argument given earlier for the 
inclusive decays, we set /if, = 2.5 GeV, and include the NNLO corrections in an analogous 
fashion to the inclusive B —> case. The explicit expressions for the B 

branching ratios can be found, for example, in Ref. |]l^ . 

The input parameters that we use in the analyses are summarized in Table Our SM 
predictions for the above discussed inclusive and exclusive branching ratios are summarized 
in Table Note that the dominant source of uncertainty comes from the form factors 
dependence. Summing the errors in quadrature we obtain: 


B{B ■ 

Ki+r) 

= (0.35 ± 0.12) X 10“® 

{5BKii = ±34%) , 

(38) 

B{B - 

^ K*e+e-) 

= (1.58 ± 0.49) X 10“® 

{5BK*ee = ±31%) , 

(39) 

B{B ^ 

K^fi+fi-) 

= (1.19 ± 0.39) X 10“® 

{6Bk*pp = ±33%) . 

(40) 


Note that the dependence of the exclusive decay branching ratios on mdmb is much milder, 
as we are using the (mc/mf,-independent) lifetime t{B^) in calculating the branching ratios 
for exclusive decays, as opposed to the inclusive decays B —> , where the semileptonic 

branching ratios are used for normalization. Since the semileptonic decay widths depend on 
me/mb, this sensitivity goes over to the inclusive decay branching ratios for B —> Xsi~^i~. 
Note also that as we have used r(R°) in calculating the branching ratios for exclusive decays, 
all the branching ratios given above are for the B^{B^) decays. The ones for the R^-decays 
can be scaled by taking into account the lifetime difference. 


5 Model independent constraints from B Xgj 

In this section we work out the 90% C.L. bounds that the measurement (P implies for 
Ay”*(2.5 GeV), where this quantity is dehned as follows: 

A*°*(2.5 GeV) = ^ 77 ( 2.5 GeV) {fiw) + A^s{2.5 GeV) Cd^{fiw) + A^^{2.5 GeV) . (41) 

It was recently pointed out in Ref. P] that the charm mass dependence of the B —>■ Xg'j 
branching ratio was underestimate in all the previous analyses. Indeed, the replacement of 
the pole mass {mc^poie/mb^poie = 0.29 ±0.02) with the MS running mass {m^^{fib)/mb,pole = 
0.22 ± 0.04) increases the branching ratio of about 11%. In order to take into account this 
additional source of uncertainty, we work out the constraints on the Wilson coefficients for 
both choices of the charm mass; we will then use the loosest bounds in the b —>■ si^i~ analysis. 
We use the numerical expression for the integrated B —>■ branching ratio as a function 
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mz 

91.1867 GeV 

as{mz) 

0.119 

mw 

80.41 GeV 

Q^e 

1/133 

'^b^pole 

4.8 GeV 

sin^ 9w 

0.23124 


(173.8 ± 5) GeV 

Gf 

1.16639 X 10-5 GeV-2 

Tbo 

1.54 ps 

\VM 

0.038 

exp 

0.104 

A 

0.225 

rric/mb 

0.29 ±0.04 

W\Vts?/W 

0.95 

Ai 

-0.2 GeV^ 

A 2 

±0.12 GeV^ 


Table 2: Input parameters and their assumed errors used in calculating the b —>■ decay 

rates. The quantities A, Ai and A 2 are, respectively, the Wolfenstein parameter and the two 
HQET parameters appearing in the heavy quark expansion. 

of Rf^sif^w) = presented in Ref. [0] (Note that, for rric/mb = 0.22, we 

had to compute the small corrections to the coefficients Bij). For the sake of definiteness we 
shall take nw = 'rnw in deriving the constraints on physics beyond the SM. We impose the 
bound (^) at 90% C.L. and include the theoretical uncertainty due to the variation of the 
scale /Xfe in the range [mb/2, 2mb]. In Fig. we present the resulting allowed regions in the 
[i? 7 (/ivi/), Rs{^w)] plane; the solid and dashed lines correspond to the rric = rric^poie and nic = 
cases respectively. According to the analysis presented in Ref. 
in Fig. to \R8{fiw)\ < 10 in order to satisfy the constraints from the decays B —>• Xgg 
and B —> (where X^ denotes any hadronic charmless hnal state). Evolving the allowed 

regions to the scale /it = 2.5 GeV and assuming that new physics only enters in we 
plot in Fig. Pd the corresponding low-scale bounds in the plane [i?7(2.5 GeV), R^i^.b GeV)], 
where i? 7 , 8 (/afe) = ^ 7 ° 8 (/^b)M 7 ^(/^fe)- The regions in Fig. Qd translate in the following allowed 
constraints: 


m , we restrict. 


r mjmb = 0.29 : A*°t(2.5 GeV) G [-0.37, -0.18] & [0.21, 0.40] , , . 

I mjmb = 0.22 : A\°\2.b GeV) G [-0.35, -0.17] k [0.25, 0.43] . ^ 

In the subsequent numerical analysis we impose the union of the above allowed ranges 

- 0.37 < A7°*’<°(2.5 GeV) < -0.17 k 0.21 < A7°*’>°(2.5 GeV) < 0.43 (43) 

calling them A 7 °*-positive and Ay^^-negative solutions. 
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B{B Xse+e-) x lO"® 


mt{GeV) 

rric/mb 

fib = 2.5 GeV 

fib = 5 GeV 

fib = 10 GeV 

168.8 

0.29 

6.30 

6.83 

7.00 

173.8 

0.29 

6.52 

7.08 

7.26 

178.8 

0.29 

6.75 

7.32 

7.52 

173.8 

0.25 

5.83 

6.30 

6.47 

173.8 

0.29 

6.52 

7.08 

7.26 

173.8 

0.33 

7.38 

8.12 

8.35 


B{B X,/i+/i-) X 10-6 


mt(GeV) 

rric/mb 

fib = 2.5 GeV 

fib = 5 GeV 

fib = 10 GeV 

168.8 

0.29 

3.70 

4.03 

4.21 

173.8 

0.29 

3.88 

4.23 

4.42 

178.8 

0.29 

4.08 

4.44 

4.64 

173.8 

0.25 

3.35 

3.70 

3.92 

173.8 

0.29 

3.88 

4.23 

4.42 

173.8 

0.33 

4.53 

4.93 

5.15 


Table 3: Dependence of the inclusive branching ratios B Xst^i [i = e, fi), in the SM on 
the scale /ib, rrit and rric/rnb- 


6 Model independent constraints from h 


In this section we compute, in the C'lo'^] plane, the bounds implied by the ex¬ 
perimental results given in Eqs. The results are summarized in Figs. ^-[T0|. In 

each hgure we focus on a different experimental bound and the two plots shown in these £g- 
ures correspond respectively to the ^^“^-negative and ^^“Dpositive solutions just discussed. 
Within each plot we then vary ^ 7 °* in the allowed range [given in Eqs. (^)]. The present 
bounds impact more strongly the decays B —>■ {X^, K*)e^e~^ for which the branching ra¬ 
tios are larger due to the smallness of the electron mass. On the other hand, the decays 
B —> Ki^£~ do not show any enhancement in the low-s region and hence they are practically 
the same for the dielectron and dimuon hnal states. Hence, the bounds for the Ke^e~ and 
cases are presented in the same plot. In Fig. 0 we combine all the bounds in a 
single plot. Note that the overall allowed region is driven by the constraints emanating from 
the decays B —>• and B . 

Some comments on the results shown in these figures are in order: 


• From the comparison of Figs. ^ and the importance of performing the analysis using 
the NNLO precision clearly emerges. In Fig. ^ we used the NLO precision ( see for 
instance in Ref. [Q). In this approximation we have to drop all the hnite corrections 
of order (that is all the terms with the superscript (1)) and the functions 007 
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B Ki+i- 

(0.35 ± 0.11 ± 0.04 ± 0.02 ± 0.0005) x 10-^ 

B K*e+e- 

(l.58 ± 0.47 ± 0.12l[(;[J^ ± 0.04) x 10"® 

B —>■ K*fi~^fi~ 

( 1.19 ± 0.36 ± 0.12^0;°® ± 0.04) x 10"® 

B Xsfi^fi~ 

(4.15 ± 0.27 ± 0.21 ± 0.62) x 10"® 

B X,e+e- 

(6.89 ± 0.37 ± 0.25 ± 0.91) x 10"® 


Table 4: SM predictions at NNLO accuracy for the various inclusive and exclusive decays 
involving the quark transition b —> For the exclusive channels the indicated errors 

correspond to variations of the form factors, fib, rnt^poie and rric/mb, respectively. For the 
inclusive channels the errors correspond, respectively, to variations of fib, mt^poie and mc/mb- 


and (Uyg; we retain the (Ug term in but drop the corresponding one in C^. The 
impact of switching on all these corrections is to lower sizably the branching ratios 
(this happens both in the full and partial NNLO scenarios previously discussed). As 
a result, the strong constraints on the new physics Wilson coefficients resulting from 
the NLO analysis are softened by the inclusion of the NNLO corrections. 

• In Fig. |T^ we identify four regions still allowed by the constraints on the branching 
ratios that present very different forward-backward asymmetries. In Fig. ^ we show 
the shape of the FB asymmetry spectrum for the SM and other three sample points. 
The distinctive features are the presence or not of a zero and global sign of the asym¬ 
metry. A rough indication of the FB asymmetry behavior is thus enough to rule out a 
large part of the parameter space that the current branching ratios can not explore. 


For the decays B Kfi^fi~ and B —> Ke^e~, a measurement is now at hand which 
we have already listed. The BELLE collaboration has combined these branching ratios, 
getting B{B = O.TSlgJi ± 0.09 x 10“® [0. In showing the constraints in 

Fig. 0 from B Ki^i , we have used this measurement to get the following bounds: 


0.38 X 10“® < B{B Ktr) < 1.2 x lO-"" at 90% C.L. 


-6 


(44) 


Concerning the upper bound, 1.2 x 10 we note that currently a discrepancy exists between 
the BELLE m and the BABAR |]^ results, with the latter reporting an upper limit B{B ^ 
< 0.5 X 10“® (at 90% C.L.) conflicting mildly with the BELLE measurements. 
However, this could just represent the vagaries of statistical fluctuations, and hopefully this 
apparent mismatch will be soon resolved with more data. Note that the branching ratio for 
B —> is bounded both from above and below, resulting in carving out an inner region 

in the (C^^(piy), <^ 10 ^) plane. 

At the end of this section we present the numerical expressions for the inclusive branching 
ratios integrated over the low-s region only where the full NNLO calculation is at hand. 
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According to the Belle analysis presented in Ref. 
follows: 


rm we choose the integration limits as 


B X,e+e" : 


0.2GeV\2 ^ /Mj/vt-0.6 GeV\^ 

) - 


mb 

\ mb / 


<s< 


\ rUb 

'Mj/^ - 0.35 GeV 
mb 


I 5 

2 


(45) 

(46) 


The integrated branching ratios have the following form: 


l3(B^XJ+e-) = 


iO 


Oi + 02 


-o vi-'y 




+04 Re Re + 05 Im Af* Im Gg^^ + Og Re A 
+07 Im + 08 Re Gg^P + 09 Im Gg^^ + oio Re Gfg’ 


■7 

rNP' 


(47) 


where the numerical value of the coefficients Oj are given in Table |] for £ = e, /r. For the 
integrated branching ratios in the SM we hnd: 


B{B Xse+e-) 
B{B ^ X,/i+/i-) 


(2.47 ± 0.40) X 10"® {6Bx,ee = ±16%) , 
(2.75 ± 0.45) X 10"® (5i3v,M/4 = ±16%) • 


(48) 

(49) 


£ 

QjI 

0-2 

Os 

04 

Os 

Qq 

OjI 

as 

Og 

Oio 

e 

1.9927 

6.9357 

0.0640 

0.5285 

0.6574 

0.2673 

-0.0586 

0.4884 

0.0095 

-0.5288 

/i 

2.3779 

6.9295 

0.0753 

0.6005 

0.7461 

0.5955 

-0.0600 

0.5828 

0.0102 

-0.6225 


Table 5: Numerical values of the coefficients Oj (evaluated at /if, = 5 GeV) for the decays 
B —>• Xs£~^£~ {£ = e, /i). Alf^ is computed at /if, = 5 GeV while G^^ at 
(Gj)()^ is scale independent). We use the full NNLO calculation which is available only in 
the low-s region. The actual ranges for the integrations are chosen according to the Belle 
analysis presented in Ref. ig. They are s G [4m^, {Mj/q, — 0.35 GeV)^] for the Xs/i+/i and 
s G [(0.2 GeV)^, {Mj/q, — 0.60 GeV)^] for the XsC+e” modes. 


7 Analysis in supersymmetry 


In this section we analyze the impact of the 6 —»• sy and b —>■ s£^£~ experimental constraints 
on several supersymmetric models. We will hrst discuss the more restricted framework of 
the minimal flavour violating MSSM, and then extend the analysis to more general models 
in which new SUSY flavour changing couplings are allowed to be non-zero for which we will 
adopt the so-called mass insertion approximation (MIA) p6t|3^ . 
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7.1 Minimal Flavour Violation 


As already known from the existing literature (see for instance Ref. [^), minimal flavour 
violating (MFV) contributions are generally too small to produce sizable effects on the Wilson 
coefficients Cg and Ciq. In the MFV scheme all the genuine new sources of flavour changing 
transitions other than the CKM matrix are switched off, and the low energy theory depends 
only on the following parameters: /r, M 2 , tan/3, Mh±, and 61 (see Appendix ^ for a 
precise definition of the various quantities). Scanning over this parameter space and taking 
into account the lower bounds on the sparticle masses (M^^ > 90 GeV, M^± > 90 GeV) 
as well as the 6 —>■ 57 constraint given in Eq. (^, we derive the ranges for the new physics 
contributions to Gg and Ciq. In order to produce bounds that can be compared with the 
model independent allowed regions plotted in Fig. we divided the surviving SUSY points 
in two sets, according to the sign of Scanning over the following parameter space 


' Mi = 90 GeV ^ 1 TeV 
Oi = -7r/2 tt/2 
tan (3 = 2.3 -7 50 
< yu = -1 TeV ^ 1 TeV 
M 2 = 0-71 TeV 
Mh± = 78.6 GeV 1 TeV 
. Mi, > 50 GeV 

we hnd that the allowed Gg and Gio ranges are: 


< 0 ^ 


e [-0.2, 0.4], 

CMFv g [_i.o,0.7] 


> 0 ^ 


e [-0.2, 0.3], 
e [-0.8, 0.5] . 


(50) 


(51) 

(52) 


We stress that the above discussion applies to any supersymmetric model with flavour uni¬ 
versal soft-breaking terms, such as minimal supergravity MSSM and gauge-mediated super- 
symmetry breaking models. Beyond-the-SM flavour violations in such models are induced 
only via renormalization group running, and are tiny. Hence, they can be described by MFV 
models discussed in this paper. 

Before finishing this subsection and starting our discussion on models with new flavour 
changing interactions, let us show in more detail the impact of b ^ s'j on MFV models. The 
scatter plot presented in Fig. ^ is obtained varying the MFV SUSY parameters according 
to the above ranges and shows the strong correlation between the values of the Wilson 
coefficients G 7 and Cg- In fact, the SUSY contributions to the magnetic and chromo- 
magnetic coefficients differ only because of colour factors and loop-functions. In Figs. 
and ^ we present the dependence of the charged Higgs and chargino contributions to Gy on 
the relevant mass parameters (that are the charged Higgs mass for the former and the lightest 
chargino and stop masses for the latter). Note that we plot the SUSY Wilson coefficients 
at the scale /Xb normalized to the SM values. In the chargino case we are able to exploit 
the 6 i and tan/3 dependence since (for non negligible values of the stop mixing angle) the 
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chargino contribution is essentially proportional to sin tan/3. Indeed, the various curves 
are very stable with respect to variations of tan (3 and 9i according to the ranges specihed in 
Eq. (1^; in order not to complicate unnecessarily the hgure we show the actual spread for the 
= 200 GeV case only (note that according to the above discussion we require \ 6 i\ > 0.01). 
In order to show the full strength of these hgures let us entertain a scenario in which Cj 
has the same sign as in the SM. In this situation large contributions to Cj are completely 
ruled out. This means that, looking at Figs. and |^, it is possible to obtain lower bounds 
on some SUSY particles. Note that Fig. ^ has very strong consequences. Assuming for 
instance = 500 GeV we see that the ratio i? 7 /(sin 6 *j-tan/ 3 ) is of order 0.2. If we 

then allow for larger values of the stop mixing angle and of tan /3, the contribution can easily 
violate the 6 —>■ sy constraint by more than one order of magnitude (e.g. for sin^j = 0.5 and 
tan/? = 50 we obtain something of order 6 that is orders of magnitude above the current 
limit). 


7.2 Gluino contributions 


Gluino contributions to Gg and Gio are governed by mass insertions in the down squark mass 
matrix. From the analysis presented in Ref. we see that the dominant diagrams involve 
the parameter {^23)ll and that large deviations from the SM are unlikely. The impact of 
{^ 23 )LR is negligible for the following two reasons. First of all, contributions to either Gg or 
Gio are obtained by —>■ sl transitions and LR insertions can therefore enter only at the 
second order in the mass insertion expansion. More importantly, the insertion (^ 23 )^/? gives 
a contribution to the coefficient Gy that is two orders of magnitude bigger than the SM one. 
The bottom line of this discussion is that (623)lr contributions to the semileptonic Wilson 
coefficients are extremely suppressed. Moreover, there are no gluino box diagrams and the 
7 -penguins are enhanced with respect to the Z-ones so that only contributions to Gg are 
non vanishing. Their explicit expression is (see Ref. for the analytical equations): 

/ 250GeV y 
\ ) 1/3 


gI’^^ = -0.93 


■( 4 ) 


23 )LL , 


(53) 


where Xgq = M|/M|, the f^^{x) loop-function is always smaller than 1/3 and can be found 
in Appendix D The situation is thus similar to the MFV case and the same conclusions 
hold. 


7.3 Chargino contributions: Extended—MFV models 


A basically different scenario arises if chargino-mediated penguin and box diagrams are 
considered. As can be inferred by Table 4 in Ref. |^, the presence of a light ^2 generally 
gives rise to large contributions to Gg and especially to Gig. In the following, we will 
concentrate on the so-called Extended MFV (EMFV) models that the two of us described 
in Ref. and that we will briefly summarize below. In these models we can fully exploit 


the impact of chargino penguins with a light t still working with a limited number of free 
parameters. 
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EMFV models are based on the heavy squarks and gluino assumption. In this framework, 
the charged Higgs and the lightest chargino and stop masses are required to be heavier 
than 100 GeV in order to satisfy the lower bounds from direct searches. The rest of the 
SUSY spectrum is assumed to be almost degenerate and heavier than 1 TeV. The lightest 
stop is almost right-handed and the stop mixing angle (which parameterizes the amount 
of the left-handed stop present in the lighter mass eigenstate) turns out to be of order 
0{mw/Mq) ~ 10%; for definiteness we will take \9^\ < vr/lO. The assumption of a heavy 
(> 1 TeV) gluino totally suppresses any possible gluino-mediated SUSY contribution to low 
energy observables. Note that even in the presence of a light gluino (i.e. Mg ~ 0(300 GeV)) 
these penguin diagrams remain suppressed due to the heavy down squarks present in the 
loop. In the MIA approach, a diagram can contribute sizably only if the inserted mass 
insertions involve the light stop. All the other diagrams require necessarily a loop with at 
least two heavy (> 1 TeV) squarks and are therefore automatically suppressed. This leaves 
us with only two unsuppressed flavour changing sources other than the GKM matrix, namely 
the mixings — 0 (denoted by and cl — 0 (denoted by We note that and 

mass insertions extracted from the up-squarks mass matrix after the diagonalization 
of the stop system and are therefore linear combinations of (^is)^/?, (^is)^^;^ and of (523 )lh5 
{ 623 )^ 1 ^, respectively. The insertions relevant to our discussion are normalized as follows: 


^u(c)i,i2 


|Um(.)I 


( 54 ) 


The phenomenological impact of ^ has been studied in Ref. |36] and its impact on the 


6 —»• 57 and b —>■ transitions is indeed negligible. Therefore, we are left with the MIA 

parameter only. Thus, the SUSY parameter space that we have to deal with is: /r, M 2 , 
tan P, Mi ^, sin %, Mh± , Mp and . 

The explicit expressions for the mass insertion contributions to the Wilson coefficients 
C 7 - CiQ are summarized in Appendix p. 

In order to explore the region in the plane (where are the sum of MFV 

and MI contributions and are explicitly defined in Appendix that is accessible to these 
models, we performed a high statistic scanning over the following FMFV parameter space 
requiring each point to survive the constraints coming from the sparticle masses lower bounds 
and h —> sy: 

‘ ^ ^ 1 TeV 

vr/lO 


Ml = 90 GeV 
9i = —vr/lO 
tan/3 = 2.3 50 

/i = -1 TeV ^ 1 TeV 
M 2 = 0-71 TeV 
Mh± = 78.6 GeV ^ 1 TeV 
Mi, > 50 GeV 

Kh = -1 -1. 


(55) 


The surviving points are shown in Fig. together with the model independent constraints. 
Note that these SUSY models can account only for a small part of the region allowed by 
the model independent analysis of current data. We stress that in our numerical analysis 
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reported here, we have used the integrated branching ratios to put constraints on the effective 
coefficients. This procedure allows multiple solutions, which can be disentangled from each 
other only with the help of both the dilepton mass spectrum and the forward-backward 
asymmetry. Only such measurements would allow us to determine the exact values and 
signs of the Wilson coefficients O 7 , Cg and Ciq. 


8 Summary 

We have presented theoretical branching ratios for the rare B decays B —> Xs£^£~ and 
B —> {K, K*)£^£~, incorporating the NNLO contributions in the former and partial NNLO 
improvements in the latter. This has allowed us to carry out a theoretical analysis of the 
radiative decays B Xgj and the mentioned semileptonic decays to the same order in 
In addition, we have included the leading power corrections in l/rub and l/rric in the inclusive 
decays. The dilepton invariant mass spectrum is calculated in the NNLO precision in the low 
dilepton invariant mass region, s < 0.25. The spectrum for s > 0.25 calculated to the same 
theoretical accuracy is not yet available. We estimate the spectrum in this range from the 
known partial NNLO, by noting that the dilepton mass spectrum in the full NNLO is close 
to the partial NNLO spectrum in the range s < 0.25 for the choice of the scale /ife = 2.5 GeV. 
Following this observation, we use the partial NNLO spectrum with this scale to estimate 
the central value of the full NNLO spectrum for s > 0.25. The branching ratios in the NNLO 
accuracy in the SM are calculated to have the values B{B —> We+e“) = (6.89 ±1.01) x 10“® 
and B{B —>■ Xs/r+/i“) = (4.15 ± 0.7) x 10“®. They are lower by typically 12% and 20%, 
respectively, compared to their NLO estimates for the central values of the input parameters, 
and are approximately a factor 2 to 4 away from their respective experimental upper limits. 
Hence, current B factory experiments will soon probe these decays at the level of the SM 
sensitivity. In view of the fact that the dilepton mass spectrum is calculated to the NNLO 
accuracy only for s < 0.25, and the long-distance effects are not expected to be dominant, 
we stress the need to measure the inclusive decays B —> Xs£^£~ in this dilepton mass range. 
In fact, as shown in this paper, such a measurement is theoretically as robust as the inclusive 
radiative decay B W 7 . 

In the second part of this paper, we have used our improved theoretical calculations to ex¬ 
tract from the current data, listed in Eqs. (I) - (0) , the allowed ranges of the effective Wilson 
coefficients - C'io(/i). In doing this, we have first determined the ranges on the Wilson 

coefficients and Cs{fi) from B —> Xgj decay, and then determined the allowed ranges 

of the coefficients Cg^{^w) and (at 90% C.L.). Since the decays B K£^£~ are now 
measured by the BELLE collaboration, they carve out an inner region in the (G^^(pvy),G(^'^) 
plane, allowed previously. Under the assumption that the SM-operator basis of the effective 
Hamiltonian is sufficient to incorporate also the beyond-the-SM physics effects, the analysis 
presented in this paper is model independent. We find that all current data are consis¬ 
tent with the SM. However, present experimental measurements allow considerable room 
for beyond-the-SM effects, which we have worked out in specific supersymmetric contexts. 
For this purpose, we have used the MFV model, and an Extended-MFV model introduced 
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in Ref. 


The resulting constraints on the supersymmetric parameters are worked out, 
in particular on the charged Higgs mass, Mh±, the lighter of the two stop masses, 
the ratio of the two Higgs vacuum expectation values, tan/3, and the MIA parameter (523)- 
With more data, expected from the leptonic and hadronic B factories, these constraints will 
become either much more stringent, pushing the supersymmetric frontier further, or else, 
more optimistically, new data may lead to impeccable evidence for new physics effects. We 
have illustrated this using the forward-backward asymmetry in R —>■ decays. 
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A Stop and chargino mass matrices 


The 2x2 stop mass matrix is given by 


M? 


M? M? \ 

tLL tLR I 

M^* M? / ’ 

tLR ^RR ' 


( 56 ) 


where 


The eigenvalues 


Mf 

^LL 

1 2 

= M? + { -sin^ 6*w) cos2/3m|-1-ruj , 

2 3 

( 67 ) 

Mf 

^RR 

2 

= M|-I-- sin^ cos 2/3 m|-t-, 

( 68 ) 

Mf 

tLR 

= mt{At - fi* cot (3) . 

( 59 ) 

are given by 


2Mf r = 

tl,t2 

+ Mfj ± ■ 

( 60 ) 


with < A37. We parametrize the mixing matrix V} so that 


\ =V'A \ +sin0£ \ ( 

U 2 / \^R j \ -sin^i cos9f j \tR j 


(61) 
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The chargino mass matrix 


M 2 mY/\/2sm. (3 
I mv(/\/2cos/3 fi 


= 


can be diagonalized by the bi-unitary transformation 

where U and V are nnitary matrices snch that M-+ are positive and M-+ < M, 

X T ^1 ' 


(62) 


(63) 


X 2 


B Wilson coefficients C 7 — Cm in EMFV models 


TO 


In this appendix we collect the explicit expressions for the Wilson coefficients C 7 - Cio in 
the mass insertion approximation. The conventions for the dehnition of the chargino mass 
matrix is snmmarized in Appendix ^ the normalization of the mass insertion is given in 
Eq. (^) in the text, and the loop fnnctions encountered below can be found in Appendix |y. 

• Contributions to the magnetic and chromo-magnetic dipole moment coefficients: 


^MI _ ^hcL 
7,8 - g 


W, 


14 


is 


Y 

M| M, ^ 

V^^Xi xMI/ 


• T/* rriiCosOi \ , 

V2smpmw ' 


sin 0iU*2 -^ (x*, ) 

mw cos p ’ 


(64) 


The contributions to the semileptonic coefficients can be divided in three classes: 
• Photon mediated penguin diagrams: 

_ PI 

n'lO — U , 


(65) 


/-iMI.'y ^ r 

^9 - 




cs 


Vt 


ts 




mediated penguin diagrams: 


^MI,Z _ 

n'lo — 


5, 


hcL 


4sin^6'iy vts iviq p 

^ ]Y-]\/[- -7 + ClT'l 2 ^- 


14, 


14 


ts 


^*2 C TA [ C /ITA* mtCosOi 


/)T/* X T/ T/* ^^2 ’ ^) 1 

-sm 0 * 14.1 4il4il42- ^ - 1 , 

/^Ad I ,Z /A * 2/1 -1 \ /^Ad I 

Cg = (4sm 0H/- l)Cio ’ 


^4 


2x 


'*2 


(67) 

( 68 ) 
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Box diagrams with an internal snentrino line: 


C 


a 


iMI.box 


Vcs 


10 

sin^ 9w 

Vts 

Ml 


X(i^^(xj, Xj, 

Xt^^Xi) , 

iMI,box 

9 — 

^MI,box 

~^10 





rut cos 6i 


v^sin/3mM.^" sin0,-lA*^ 


(69) 

(70) 


The branching ratios for the various decays are obtained from Eqs. (0-® by means of 
the following replacement: 


V I /^IVll 

'^7,8 '-"7,8 “T '-"7,8 5 

9,10 


yMFV I j-yMI 


CNP 


^MFV I I ^MI,Z . ^Ml^box 

^ 9,10 “T <^ 9,10 “T ^ 9,10 “T ^ 9,10 


where the expressions for can be found in Ref. . 

C Loop functions 

The various loop functions introduced in Appendix are listed below. 

—7 -I- 12a: -|- 3a:^ — 8a:^ -|- 6a:(—2 -|- 3x) logx 


h{x) = 
f2{x) = 
fsix) = 
hix) = 
= 


6(x- 1)4 

5 — 12x -I- 7x^ — 2x(—2 -I- 3x) logx 
2(x- 1)3 

2 -|- 3x — 6x^ -|- x3 + 6x log X 
6(x-l)4 ’ 

— 1 -1- x^ — 2x logx 
2(x- 1)3 ’ 

/j(l/x) — fiiy/x) 

^ 1,2,3,4) . 


fsix) 


L^\x,y) = 


j{x) 

k{x) 


X logx 

X — 1 

.2 


x(l - y) 

52 — 153x -I- 144x^ — 43x3 -|- 6(6 — 9x -|- 2x3) logx 
6(x- 1)4 

2 — 9x -I- 18x^ — 11x3 -f 6x3 jQg^ 

(x — 1)4 
fi{x) - fi{x/y) 


1-2/ 


(7 = 7,8) 


.. ^ j{x)-j{y) j{x,z) - j{y,z) 

J{x, y) = -, 3 [x, y,z) = - 


X^ log X 
X — 1 


k{x,y) = 


x-y 

k{x) - k(y) 
x-y 
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Hx,y,z) = 


x-y 

k(x,z) - k(y,z) 
x-y 


(71) 

(72) 


(73) 

(74) 

(75) 

(76) 

(77) 

(78) 

(79) 

(80) 

(81) 

(82) 
























d 2 {x,y,z,t) 

d^\x,y,z,t) 


X logo; 


{x — y){x — z){x — t) 
d2{x,y, l,t) - d2{x,y,z,t) 


+ {x y) + {x ^ z) + {x ^ t) 


, (83) 
(84) 
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Figure 1: Relative size R{s) of the eombined l/rrih and l/rric power corrections as defined 
in Eq. in the decay rate in B ^ decays as a function of the dilepton invariant 

mass in the SM (solid) and for C-j = (dashed). 




Figure 2: Partial (dashed lines) vs full (solid lines) NNLO computation of the branching 
ratio B —>• In the left plot (s G [0,0.05]j the lowest curves are for /i = 10 GeV 

and the uppermost ones for p, = 2.5 GeV. In the right plot the p dependence is reversed: 
the uppermost curves correspond to p = 10 GeV and the lowest ones to p = 2.5 GeV. The 
right-hand plot also holds for the decay B —>• Xgp'^p~. 
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Figure 3: 90% C.L. bounds in the plane following from the world average 

B —>■ Xsj branching ratio for /i = mw (left-hand plot) and ja = 2.5 GeV (right-hand plot). 
Theoretical uncertainties are taken into account. The solid and dashed lines correspond to 
the uic = nic^poie and rUc = cases respectively. The scatter points correspond to the 

expectation in MFV models (the ranges of the SUSY parameters are specified in the text). 




Figure 4: NLO Case. Constraints in the [Cg^{fiw), plane that come from the BELLE 
90% C.L. upper limit B{B —>■ Xse'^e~) < 10.1 x 10“®. Theoretical uncertainties are taken 
into account. The plots correspond to the A\°^{2.5 GeV) < and CeF) > 0 case, 

respectively. In each plot the outer contour corresponds to the smaller \ A\°^\ value. The dot 
in plot on the left is the SM point. 
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Figure 5: NNLO Case. Constraints in the [Cg^{fj,w),CiQ^] plane that come from the 
BELLE 90% C.L. upper limit B{B —> X^e+e”) < 10.1 x 10“®. See Eig. ^ for further details. 




Figure 6: NNLO Case. Constraints in the [Cg^{p,w),CiQ^] plane that come from the 
BELLE 90% C.L. upper limit B{B < 19.1 x 10“®. See Eig. for further details. 
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Figure 7; NNLO Case. Constraints in the [Cg^{fiw),CiQ^] plane that come from the 
90%CL BELLE constraints 0.38 x 10“® < B{B —>■ Ki^i~) < 1.2 x 10“® (Eq. [^. Note that 
only the annular regions between the two circles is allowed. See Eig. ^ for further details. 




Figure 8: NNLO Case. Constraints in the [Cq^{p,w),Ciq^] plane that come from the 
90%CL BELLE constraint B{B —> iF*yU+/r“) < 3.0 x 10“®. See Eig. for further details. 
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Figure 9; NNLO Case. Constraints in the plane that come from the 

90%CL BELLE constraint B{B —>• K*e^e~) < 5.1 x 10“®. See Eig. ^ for further details. 
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Figure 10: NNLO Case. Superposition of all the constraints. The plots correspond to the 
GeV) < 0 and GeV) > 0 case, respectively. The points are obtained by 

means of a scanning over the EMEV parameter space and reguiring the experimental bound 
from B —>• Xs 7 to be satisfied. 
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Figure 11: Differential Forward-Backward asymmetry for the decay B Xgt^l . The four 
curves correspond to the points indicated in Fig. [7^ . 



Figure 12: Dependence of Rj^{jab) = C^^{pLb)/Cj^{pbb) on the mass of the charged Higgs. 
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Figure 13: Dependence of RJ{f^b) = Cj{nb)/C^^{inb) on the mass of the lightest stop in MFV 
models. The chargino contribution is essentially proportional to sin6*ftan/9 for not too small 
sin6'j. For the set of curves 2 we show the variation due to several choices of 9^ and tan p. The 
thick bunch of lines is obtained for (sin6^^, tan/9) = (0.025,40), (0.05,20), (0.1,10), (0.2,5), 
(.5,2) for which sin 6*^tan= 1. The two thin curves correspond to (sin 6*j, tan/9) = (0.1,20) 
and (.1,40), with the product sin 6*^tan having a value 2 and 4, respectively. 
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